ON THE INTEGRAL OF THE PRODUCT OF FOUR 
AND MORE BERNOULLI POLYNOMIALS 



DAEYEOUL KIM, SU HU, AND MIN-SOO KIM 

Abstract. L.J. Mordell in [Integral formulae of arithmetical char- 
acter, J. London Math. Soc. 33 (1958) 371-375] provided the for- 
mula for the integral of the product of two Bernoulli polynomials, 
he also remarked: "The integrals containing the product of more 
than two Bernoulli polynomials do not appear to lead to simple 
results." Let r be any positive integer. In this paper, we provide 
the explicit formulas for the integral of the product of r Bernoulli 
polynomials. Many authors' results in this direction, including 
Norlund, Mordell, Carlitz, Agoh and Dilcher are special cases of 
the formulas given in this paper. 



1. Introduction 

Nearly 90 years ago, Norlund gave the following formula for the 
integral of the product of two Bernoulli polynomials in his book "Vor- 
lesungen uber Differenzenrechnung" : 

(1.1) / B^B^dz = (-l) k - 1 ——B k+l , k + l>2 

Jo {k + iy. 

(see [HI p. 31]). 

Nikolas and Mordell provided two different proofs of (11. ip in [8] 
and [7] respectively. Mordell remarked: " The integrals containing the 
product of more than two Bernoulli polynomials do not appear to lead 
to simple results." (See [H p. 375]). But Carlitz [3] obtained an explicit 
formula for the integral of the product of three and four Bernoulli poly- 
nomials which extends (II. ip . Later, Wilson [10] generalized Carlitz's 
result on the integral of the product of three Bernoulli polynomials by 
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evaluating the following integral 



1.2) / B k (az)B l (bz)B m (cz)dz, 



o 



where Bk(x) is the periodic extension of Bk(x) on [0, 1) and a, b, c are 
pairwise coprime integers. Carlitz's result became special case when 
a = b = c = 1. Similar integral evaluations have also been used by 
Espinosa and Moll [5] during their study on Tornheim's double sums. 

Recently, Wilson's result has been generalized by Agoh and Dilcher 
[2]. In fact, they proved the following result: 



Proposition 1.1 ([21 Proposition 3]). For k,l,m>0 

k+l a 

£(-D"£ 

where 



Ul\m\ ^ ^ \i ) (k - a + - i)\{m + a + 1)! 



X 



h,i,m = / B k (z)Bi(z)B m (z)dz, 
Jo 

C k ,i, m = B k (x)Bi{x)B m {x) - B k BiB m . 

Let r be any positive integer. In this paper, we provide the explicit 
formulas for the integral of the product of r Bernoulli polynomials 
which lead all the previous mentioned results for the integral of the 
products of Bernoulli polynomials become special cases (see Remark [T75l 
below) . 

Before statement our main results, we introduce the definition and 
some basic properties on multinomial coefficients. Since the Binomi- 
als coefficients have been used by Agoh and Dilcher in [2] to get their 
formulas of the product of three Bernoulli polynomials, in order to get 
more general ones, we need some information on multinomial coeffi- 
cients. They are defined to be the number: 



l • • • 4 kj> J k\ ! * * * kf\ 
where k\ + • • • + k r — fx and k±, . . . , k r > 0. When r = 2, we obtain the 
binomial coefficients 



(1.4) 



ki,k 2 ) \kufi-ki) \ki) \k 2 



The multinomial coefficient can be extended as follows 

w \ |0 if minjfcj} < and max{/cj} > /i 

^ j ^ j *' ' l I I < < ; ' 1 < 



l<i<r l<i<r 

fci, . . . , k r J |^ | otherwise. 
Lemma 1.2. Multinomial coefficients have the recurrence relation 

k\, . . . , k r J \k\ 1, . . . ,k r J \k\, . . . , k r 1 

and symmetry 

V \ = ( V \ 

K k\, . . . , k r J \k,j^, . . . , k a (r)J 

where a is a permutation of (1, ... ,r). 

The following two are the main results of this paper. 
Theorem 1.3. Let 

h u ...,kr = / B kl (z) ■ ■ ■ B kr (z)dz, 
Jo 

Ck u ...,k r = B kl (x) ■ ■ ■ B kr (x) — B kl ■ ■ ■ B kr , 
(1.6) ? 1 j 

J fci,...,fc r ^ ! • • • jfc J kl,---,kr 1 

1 

jfc 1 !...fc r ! Ufcl, "" fcr ' 
For any /x > 1 and fci, . . . , fc r > 0, we /love 



Cfci,...,fc r — "r - : i—;C k 



(1.7) 



? < * = E<-i)- E % ,) 

a=0 ii+--H r -i=a 
0<«i,...,i r -i<a 

^ Cfci— ii,...,fe r _i— i r _i,A; r +a+l 

+ E (■ ^ ■ ) 



iiH h«r-i=M 

0<ii,...,i r -i<fi 



X /, 



Letting 



'1.8) h 1 ,...,k r = when min < 0. 

Ki<r 



and 



fi — ki-\ h A; r _i + 1 
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in Theorem 11.31 we obtain the following result. 
Theorem 1.4. For any ki, . . . , k r > 0, we have 

fclH hfe r _i , v 

? - -= E £ U % J 

a=0 ti+--Hr-l=a n- 1 '" - ' ' l ' 

0<ii,...,i r — i<a 

Remark 1.5. In Theorem 11.41 when r = 2 and x = 1, we obtain the 
results of Norlund [9] and Mordell [7], when r = 3 and x = 1, we obtain 
the result of Carlitz [3J, when r = 3, we obtain the result of Agoh and 
Dilcher [Tj. 

Using Theorem 11.31 we immediately get some concrete examples for 
the integral of the product of four Bernoulli polynomials. 

Example 1.6. (1) I liljlfl = \B\B 2 + \B 4 - \B 2 . 

(2) h, 1,1,3 = + ^Bq - |S 4 . 

(3) A, 1,1,5 = 2 BfB 6 + 56^8 - 12^6- 

(4) h, 2,2,1 = 2-^2 + 2^2-84 + -B1-B4 + g-B6 — e-B4- 

(5) 4,3,3,4 = — YqB^Bs — -^B 2 B W — Y^BfB 10 — 495-812 + -^B w . 

From Theorem II .41 and the definition of the integral Ik u ... t kri we have 
the following result. 



Corollary 1.7. Let ^(1) 6e t/ie right-hand side of Theorem I.4 

and let a G S r , where S r is the symmetric group of degree r. Then for 
all ki, . . . , k r > 0, 

T kl ,...,k r {x) = T a ( kl ) t _ )a ( kr )(x). 

The following parts of this paper are organized as follows. 

In section 2, we recall the definitions and some properties for the 
Bernoulli polynomials, in particular the well-known fact for the integral 
of one Bernoulli polynomial. They will be used in the proof of the main 
result of this paper. In sections 3 and 4, we prove Theorems 11.31 and 
11.41 respectively. In the last section, by using our formulas, we derive 
several examples including some formulas by previous authors. 



2. Bernoulli polynomials 

The Bernoulli polynomials B k (x) are denned by the exponential gen- 
erating function 

fp xt J° fk 

k=0 

They satisfy the following well-known identities: 

(2.2) B k (x + 1) - B k {x) = kx k ~\ k > 1, 

(2.3) 4~B k (x) = kB k ^(x), k > 1, 
ax 

rx+l 

(2.4) / B k (z)dz = x k , k> 0, 

J X 

(2.5) jy t (z)d,= B ^i-^ x \ k>o, 

(2.6) B k (\ - x) = (-l) k B k (x), k > 0. 

Let B k = B k (0) be the Bernoulli numbers. From the above definition 
of B k (x), we conclude that all the B k are rational numbers. It is well- 
known that B 2 k+i = for k > 1, and B k (x) are alternatively positive 
and negative for even /c. Letting x = in (12. 6p . we get the following 
special value of the Bernoulli polynomials: 



(2.7) B k {l) 



B k for jfe ^ 1 
Si for = 1 



(see 0). From (J2l} and (l2"31) . we have 



.•_n „-_n •/ i 

(2.* 



i=0 i=0 

B k (z)dz 

o 

_ B k+1 (n + 1) - gfc+i 

fc+ 1 ' ~ ' 

This is the formula for the integral of one Bernoulli polynomial. 
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3. Proof of Theorem 11.31 

Using integration by parts, along with (11 ,6p and (12. 3p . we obtain 

If d 
h u ...,k r = ^— J B kl (z) ■ ■ ■ B kr _ 1 {z)—B kr ^ 1 {z)dz 

= [B kl {z) ■ ■ ■ B kr _ x {z)B kr+1 (z)]l 

(3-1) ~KTlJ ^- z (B kl (z)---B kr _ 1 (z))B kr+1 (z)dz 

1 ° 

< ^ y fcl,...,fc r _l,fc r + l 



k r + 1 
~ k 



Tl [ Tz {Bkl{z) ' ' ' B ^))B kr+ i(z)dz. 



From (12.31) we have 



±{B kl {z) ■ ■ ■ B^z)) 
( 3 - 2 ) =k 1 B kl _ 1 (z)---B kr _ 1 (z) 

H h k r _ x B kl {z) ■ ••B kr _ 1 _ 1 (z). 

From above and (11.61) we have 

i r d 



k r + 1 ,/n dz 

h 



(3.3) /,V 



(B kl (z) ■ ■ ■ B^z^B^+^dz 

- / B kl _ x (z) ■ ■ ■ B kr ^(z)B kr+1 (z)dz 
1 Jo 



+ ■■■ + ^7 I B kl (z)..-B kr _ 1 -. 1 (z)B kr+1 {z)dz 

&r + J- JO 

J fcl-l,...,fc r _i,fc r + l T~ T " ~ _ r-'A:i,...,fc r ._i-l,A; r + l; 



k r + l — fcr + 1 

and substituting ( 13. 3f) to (13. lj) we get 
(3.4) i 

Ik\,...,k r | rCA,'i,...,fe r _i,fc r +l 

rC r . "T" 1 

~~ 1 ~T(^1^1-1 r-Xi* r+1 + ' ' ' + ^r-l-^fci,...,fc r _i-l,fc r + l)- 

fc r + 1 

From (11.61) and (13. 4p . we have 
(3.5) 

Iki,...,k r Cfci,...,fc r — i,fe r +l \J-ki — l,...,fe r -i,fc r +l "I - ' ' ' Ik\,...,k r -\ — l,fc r +l) 
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Now we prove the case when /i > 1 in Theorem 11.31 by induction. For 
\i = 1 it just reduces to (13. 5p . Now we assume that Theorem 11.31 holds 
for some fj,>l. Write 



s = E (■■)'■ 

/ — ' Wi, . . . ,z r _i I 



(3.6) 

HH hV-l=A» 

0<i\,...,i T —\<fx 



Eq. (11.71) becomes 

'* - = B-ir E (,,..°,J 

a=0 ii+-+i r _i=a v " ) r 1/ 

(3 7) 0<ii,...,i r _i<a 

X Cfci— ii,...,fe r _i— i r _i,fc r +a+l 

+ (-lfS. 

From (13.51) we have 

S= ^ ^ (. . ) Cfe 1 -ii,...,fc r _i-i r _i,fe r .+/i+l 



iiH Nr-i=A t 

0<ii,...,i r _i</i 



(3.8) u+-+v-i=i" 

0<ii,...,i r _i</i 



E ■ ■ * 



ki — (ii + l),...,k r -i—i r -i,k r +fJ.+l 



E ■ * 

tlH Mr-1=M 

0<ti,...,i r _i</i 



fei— ji,...,fc r _i-(i r _i + l),fc r +^+l- 



Let = ij + 1 in (I3.8P for j = 1, . . . , r — 1. From (11.51) . we obtain 

(., / • )=o if<i = o, 

- 1, • • .,« r -l/ 



(3.9) 







if = 0, 



s 
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and 



A* 

i'i 1, . . . , i r -i 



= if i 2 = A* + 1, • • • , V-i = A* + 1> 

V*i ■•■}•••) V— 1/ 
(3.10) 

= if i x = a* + 1, . . . , z r _ 2 = A* + 1; 







then we have 



iiH Kr-i=£« 

0<ii,...,i r _i</i 



E ■ ■ * 

^ V «i,..., v-i 7 



fcl— (ii+l),...,fc r _i— i r _i,fe r +M+l 



^ W'x - 1, . . .,l r -lj 



(3.11) 



i^H h«r-l=M+l 

0<i' 1 ,...,i r -i<At+l 



E (■ 



iiH hv- 1=/-« 

0<«i,...,i r _i</i 



fci-ii,...,fc r _i — (i r _i+l),fc r +/^+l 



E 



\-i' r _ l =H+l 

0<ii,...,i;_ 1 </*4-l 



A* 

«1, • • .,*r-i 



'fci-ii,...,A; r _i-iJ,_ 1 ,fc r +/.t+l • 



From (J3ISD and fl3TTTD we get 



iiH hir-i=M 

0<ii,...,i r _i</i 



^ 5 v ( • ■ ) Cfcl— n,...,fc r -l— ir-l,fer+W+l 

^ . . . ,l r -lj 



(Qio) »iH h«r-l=A»+l 

1 ' o<i;,...,v_i<M+i 



^ ] ( ■/ -i • ) ^fci-i' 1 ,...,fc r _i-i r _i,fc r +^+l 

*X 1 1 ■ ■ ■ i *r— 1/ 



(• ■/ i ) -^fel-il,...,fcr-l-ir-l> fe r-+M+l- 



ixH h«J._ 1 =A»+l 

o<n,...,i;_ 1 </Lt+i 



Let ij = i'j for j = 1 , . . . , r — 1 in (I3.12p . From Lemma 11.21 (I3.12p 

becomes 
(3.13) 

S = y I . . 1 C^—f, kr-i— tr-i.fcr+u+l 



llH h«r-l=A t 

0<ii ,...,i r _i</i 



E 



iiH h«r-l=A t +l 

0<ii ,...,i r _i</i+l 



■ / ■ + "' + - " , 

*1 - 1, • • • , l r -Xj \H, ■ ■ • , V-l ~ 1 



^ ^ki — ii,...,k, — i— i, — i,fer+^+l 
/' 



iiH Nr-1=A» 

0<u,...,i r -l<A» 



/ I . • I ^fci— ii,...,fe r _i— i r _i,fc r +/i+l 

^ V zi, . . . ,i r -ij 



E (■ " +1 )'~ 

UH Mr-1=^+1 

0<ii,...,i r -i<M+l 



fci— ii,...,fe r _X— i t — i,fc r +^+l • 



Substituting (13 . 131) to (13. 7ft . we get Theorem 11.31 when fi + 1 which 
completes our proof. 

4. Proof of Theorem 11.41 
Letting 4 1 -i 1 ,...,fc r _ 1 -v_i,fc r + /i 7^ in TheoremOl from (TO]) , we have 
(4.1) ki-ii>0,..., k r _i - i r _ x > 0. 

Thus 

(kx-ii)A h (*v-i - V-i) 

= fci H + fc r -i — (ix H h v_i) 

(4.2) 

= «i H h A; r -i - // 

> 0. 

Let fj, — kx H h fcr-i + 1 in (14. 2p . We have fci H h fcr-i - // < 0. 

This contradicts with (14. ip . Therefore when [i = kx + • • • + k r -x + 1, 
one of 

kx ~ i\i ■ ■ ■ ■> k r —i ~~ V— l 
will be always negative. Thus from (11.81) . any terms 

^fei — i\,...,k T —\— i, — i,feiH hfci — i+Av+1 



10 DAEYEOUL KIM, SU HU, AND MIN-SOO KIM 

vanishes. The result then follows from Theorem 11.31 



5. Some further consequences 



Letting r = 2 in Theorem 11.41 we have 

fci 

-ffei,fe 2 = l) a Cfci-a,fc2+q+l- 
a=0 

This is equivalent to the following proposition: 
Proposition 5.1 ([1, Proposition 1]). 

I B " (s)B " w<fa = (*, + *!, + i)i S ( ) V J 

x (^Jfei-aW^fca+a+lW — B kl _ a B k2+a+1 



Letting r = 3 in Theorem 11.41 we have: 

Proposition 5.2 ([U Proposition 3]). 

fci+fc 2 



^fci,fc 2 ,fc 3 — /_, (~ I)" /, I • ■ ) Cfei-ii,fe 2 -i2,fc3+a+l 

a=0 ii+i 2 =a 
0<ii,i 2 <a 

fci+fc 2 a /_\ ^ 

^fcl — a+i,fc 2 — i,A:3+a+l 



X>ir£ 

a=0 i=0 



z/ (ki - a + i)\(k 2 - i)\(k 3 + a + 1)! 



Letting r = 4 in Theorem 11.41 we have: 



Proposition 5.3. 

fcl+fc 2 +£)3 

, .... / " W 

y fei— ii,fc 2 — 12, k3— «3,A;4+a+l ■ 



fcl+fc 2 +fc3 y v 

h M3 M= E (- x ) a E LtJ^ 



a=0 «i+« 2 +j3=a 
0<ii,i 2 ,i3<a 



Let 



Write 



hiMMfoiX) = \ B kl (z)B k2 (z)B k . A (z)B ki (z)dz. 
Jo 



~ 1 
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Using fll.6p . (12. 6 p and Proposition 15.31 we reobtain the following known 
evaluation of the integral from to 1: 

E E 



(i) 



l k 1 ,k 2 ,k s ,k4 



hlk^M 



a=0 



ii+«2+*3=a 
0<ii,i2,i3<o 



a 



X ^_l^i+k 2 +k3+k 4 +l _ j\ 



SO 



l k 1: k2,k 3 ,k4 i \a+l d 

a=0 

E 



^1+^2+^3 



(5.1) 



+i 



x 



H +12+13=0 

0<ii,i2,i3<a 



a \ ~ ~ ~ 

• ) ^ki—h Bk 2 -t2 ^k 3 — i 3 
«1 9 «2,«3/ 



if fci + & 2 + fc 3 + fc 4 = (mod 2), and 







if ki + &; 2 + k 3 + A; 4 ^ (mod 2). Therefore, without loss of generality 
we assume that 

(5.3) fci + A; 2 + A;3 + £; 4 = (mod 2) and & 4 + a + l = (mod 2) 

since B^+a+i = 0if£; 4 + a + 1^0 (mod 2). By ( 15. 3ft . we see that 
(5.4) 

fci + fc 2 + fc 3 + ^4 + 1 = (h - ii) + (k 2 - i 2 ) + (A* - i 3 ) + {h + a + 1) 

= 1 (mod 2). 
Thus from (15 .ip and (15.41) we get the following table: 





h - ii 




h - i 3 


k 4 + a + 1 


(A) 


odd 


even 


even 


even 


(B) 


even 


odd 


even 


even 


(C) 


even 


even 


odd 


even 


(D) 


odd 


odd 


odd 


even 



From now on, we will use the fact that £? 2 i+i = for i > 1. Thus, if 
kj — ij are odd for j = 1, 2, 3, then fcj — ij must be 1, so we only need 
to consider the following four cases: 
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Case (A). We put i\ = k\ — 1. Then by (15. 3p and a = i\ + i 2 + 23, 
we have 

k 2 -i 2 = k 3 -i 3 (mod 2). 
Case (B). We put i 2 = k 2 — 1. Similarly, we have 

ki — i\ = k 3 — i 3 (mod 2). 
Case (C). We put i 3 = k 3 — 1. Similarly, we have 

ki — i\ = k 2 — z 2 (mod 2). 

Case (D). We have 

h = k x - 1, i 2 = k 2 - 1, h = h - 1 

since i? 2 i+i = for i > 1. 

Using ([51]), O, (A), (B), (C) and (D), we obtain the following 
theorem. 

Theorem 5.4. (1) 7/ k x + fc 2 + A; 3 + k A = (mod 2), tfien 

( — l) a 7?fc 4 +a+l 



(1) 

fciljyjfesliy 



a=0 

\ a-fci+1 / ; i i\ 

a \ < — , / a — k:-, -+- I \ 

Bbo-iB, 



1 7 i=0 

x a— fo+l / , -, 

a \ \ -\ / a — K2 + 1 



k$— i-L>k\+k2+i— a— 1 



b, ■) e r ; 



i=Q 

a \ / a — K3 + 1 



fci — i^k2-\-k%+i— a— 1 



i=0 



k\ — i-^k2+k^+i— a— 1 



+ i(_ 1 )fci+fc2+fe Ai + fc 2 + ^3 - 3\ (k 2 + A; 3 - 2 



X 7?/ Cl+ fe 2+ fc 3+ fc 4 _2- 

(2) If k x + A; 2 + h + k A ^ (mod 2), tfiera 

(!) n 

Remark 5.5. Our formula in Theorem 15.41 is equivalent to (5) in Carlitz 
Pp. 361]. 
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Letting k\ = and ki = k, k% = I, = m in Theorem 15.41 we obtain 
the following formula for the integral of the product of three Bernoulli 
numbers: 

Corollary 5.6. For all k,l,m> 1, 

r (-ir +i kw. m \E k a =o [G_\) + (A)] wA+z-a-i 

4,z,m(l) = < ifk + l + m = (mod 2) 

[o ifk + l + m=£0 (mod 2). 
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